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ABSTRACT 


As part of his management planning and  ccntrol 
function, the U. S. Coast Guard's Pacific Area Commander 
schedules the operational missions for all Hign Endurance 
Cutters in the Pacific Area. To provide a powerful 
Management tool to assist this scheduling process, an 
analytic model for this large scale problem has been 
developed and implemented. It contains nission requirements, 
restricted Sequencing of missions, Ships! physical 
limitations and crews! morale-related considerations. The 
modeling approach is based on the Geofírion-Graves model for 
parallel production lines with significant changeover costs. 
The implementation solves a large (860 гом) 
Koopmans-3eckmann ixed charge Quadratic Assignment model 
using a new method with an advanced, feasible starting 
Solution provided by an inbedded network (with 1,720 nodes 
and 739,500 arcs). Many linear programming problems (200 
row, 450 variable) are then solved with a linear programming 
subroutine of advanced design. The resulting model and these 
implementation techniques produce excellent quality working 
schedules with very reasonable execution time and memory 


requirements. Alternative solutions are easily generated. 
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I. INTRODUCTION 


ee Seg ee Dom emp EE сш озшш шы ш 


А. BACKGROUND 


The United States Coast Guard is responsible for many 
areas of national concern in the maritime regions of this 
country. These responsibilities, called "missions" ру the 
Coast Guard, include search and rescue, law enforcement, 
fishery and custom regulation, small craft and commercial 
vessel safety, icebreaking, International Ice Patrol and 
collection of oceanographic data. To coordinate these 
responsibilities, two major commands have been established -~ 
ве for the Atlantic, Gulf of Mexico area, and one for the 
Pacific area. Each area is further subdivided into several 
HN UScts. The District Commanders handle all operational 
and administrative matters arising in their districts. The 
Area Commander, however, assumeS operational control for all 
Missions that are "multi-district" in geographical region, 
SCCpe, Or resources. Some of these misSions are: law 
enforcement and fishery patrols in the Alaska region, 


oceanographic data collection, and readiness training. 


To fulfill bis responsibilities, the area Commander must 
supervise the allocation of resources for these missions. 
Periodically, the Pacific Area Commander's staff prepares 
and formally issues the operational schedules for the High 
and Medium Endurance Cutters in the Pacific Area. This 
published schedule also satisfies tne maintenance, military 
Beadiness, and training requirements. It serves as the 


Operational orders to the Commanding Officers of tne ships, 
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specifying where, for what, and when the ships are needed. 
It informs the crews when they can expect to be home with 
their families and wnen and for how long they NM be away. 
The Coast Guard's greatest asset is its people; the schedule 
has a large impact upon their lives. Thus, the published 
schedule is cf prime concern to all - the Area Commander, 
District Commander, their staffs, ship's Commanding Officer, 


and the manning personnel and families. 


BemeeGURRENT SCHEDULING METHOD 


In any organization, the scheduling of resources is part 
Se tne Managerial planning and control function. The areas 
of responsibility of the Coast Guard, by their nature, 
generate reguirements that Can never be completely 


Satisfied. Law enforcement and protection of lives and 


property ars examples of missions that can never pe 
over-performed. The Coast Guard possesses a finite number 
of resources - ships and personnel. These resources have 


physical limitations and needs that limit their utilization 
to fulfill the mission requirements. The managerial planning 
and control function of the Area Commander is to allocate 
his resources, striking a balance between the requirements 


aaa echese constraints. 


To furnish the schedulers with basic information 
concerning what missions are required, what resources are 
available, and what limitations exist, the Area Commander 
has estaolishedä a set of guidelines. These guidelines 
compile Coast Guard-wide, Area, and District Осу 
decisions and plans by specifying patrol standards, lengths 
and frequencies for each mission. Resource limitations are 
expressed as annual maintenance and training requirements. 


The needs of the manning personnel are reflected in workload 


MZ 





balancing, Maximum cruise lengths, minimal times to be spent 
in homeport, and other items that affect the crew's 


morale". 


These morale-related limitations result in goals that 
the scheduler attempts to satisfy. The evaluation and 
comparison of two different schedules that satisfy the 
mission requirements and morale-related goals in different 
balances can not be performed using strictly objective 


criteria. Judgement and individual preferences are involved. 


Using the guidelines, two schedulers spend about one 
week manipulating a large magnetic Gantt chart and 
calculating summary totals for each candidate schedule (e.g. 
total time away from homeport, total travel time for each 
Ship). Then, additional time is spent with the resultant 
Schedule in negotiations between the decision makers, the 
Schedulers, and the Coumanding Officers of the ships. The 
desired final schedule is one that satisfies the mission 
requirements with consideration and respect for the 


morale-related goals. 


Each January and July a formal schedule is published 
covering the next 18 month period. A tentative scaedule for 
tne subsequent six month period (months 19 to 24) is also 
prepared, but distributed informally. The final published 
schedule provides one-day time zeselut:on, Specifying 
required  on-scene arrival and departure dates and estimated 


ncmeport departure and return dates. 


COAL OF TEKIS STUDY 


The inmediate goal of this exploratory research is the 


development of an analytic scheduling model based on 
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mathematical programming techniques. This analytic model, 
incorporating the most important of the mission requirements 
and morale-related goals, will generate potential working 
Schedules using a reasonable amount of computer resources. 
These potential schedules will present the decision makers 
and schedulers with more choices СОЕ objective and 
subjective evaluation than the current manual method. This 
managerial tool will aid strategic and tactical decision 
processes, and possibly improve them by relieving a heavy 


clerical worklcad. 


mrstOraically, scheduling problems (as an application of 
mathmatical programming) have either been oversimplified 
with a resultant loss of effectiveness, or so detailed that 
Semputational difficulties prevent economical solution. 
Certain recent advances in optimization capabiliries, 
however, have made possible this exploratory research to 
apply mathematical programming models to the Coast Guard 


Scheduling problem. 
The chapters of this thesis follow a systems anaiysis 
approach: 
1. Formulation (identify problem and objectives) , 


2. Analysis (examine eleuents, interrelatıonships, 


variables, and constraints), 


3. Synthesis (classify the problem and study alternative 


metncds of solution), 
4. Selection of Method and Detailed Modeling, 
5. Implementation (the test of all that precedes), and 


б. Evaluation and reformulation. 
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В ПАГ ОГНЕ AND SCHEDULING GUIDELINES 


Oo -— a a ete Í -E e —S eee eee eee тығы -— сшщ Ep oov. сш A шшщщ | | сш эшш A шша ш хш шшш _ 


EN  u:DULE'S MEASURE OF EFFECTIVENZSS 


Шеке typical production or job shop scheduling, the 
Coast Guard problem does not possess a single measure of 
effectiveness on which objective evaluation of candidate 
schedules can easily be based. There are some objective 
evaluation eriter3a, However, the subjective criteria 
derived from the morale-related considerations play tne 
dominant role in evaluation.  Calendar-related events are 
also considerations (e.g. fishing seasons, major ncliday 
periods, and weekends). Individually, each requirement or 
morale-related goal is reasonable. Unfortunately, there is 
no schedule that can simultaneously meet all the desired 
goals. Since any scheduie will have violations of the 
guidelines, it iS necessary to assess the trade-offs among 
the conflicting requirements, and to balance the various 
goals. Managerial judgement and personal preferences (of 
the Area and District Commanders, their staffs, and each 
ship's Commanding Orne om) ultimately determine the 


properties and structure of tne final schedule. 


BASIC ELEMENTS 


There are three basic elements in the schedule - ships, 
time, апа missions. There are presently Seven High 


Endurance Cutters (HEC) and five Medium Endurance Cutters 
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(MEC) in the Paciric Area that are scheduled by the Area 
Commander. cae” Carra ves 1 June 1977. The 
"Hamilton class" of HEC's, because of their different 
capabilities, are used for different types of missions than 
the single older HEC in the Area. Scheduling interactions 
between these ship types and all other ships are minimal. 
Thus, the scheduling process can be divided into two 
separate cases - one for the "Hamilton" class HEC's and one 
for the other, smaller cutters. то Fl init scope; and 
emphasize basic issues, this study concentrates cn the 
Ни 1129 ргосезз гог just the "Hamilton" diga Endurance 
Cutters (HEC's). Two of these HEC's are based (homeported) 
momaonOolulu, Hawaii; two in San Francisco, California; and 


the rest in Seattle, Washington. 


As mentioned above, each published schedule covers ап 
eighteen month period. Months 19 to 24 are tentatively 
Scheduled but not fornally distributed. Between successive 
SChedules, there is a desire to minimize the number of major 
changes. Unless major changes in mission requirements have 
occurred cr unforeseen events have taken place, the first 
Six months of a new Schedule will be in close agreement with 
the same period for the previous schedule. As the time 
frame Of the schedule progresses, more changes and thus 
greater deviations occur. The informal 19 to 24 month 
segment is the initial projection for tnis calendar period; 


Bra, therefore, highly speculative. 


The missions performed by the "Hamilton" class High 


Bus ucc Cutters cover most of the Coast Guard's areas oz 


responsibility - fishery and custom regulation, law 
enforcement, oceanography, search and rescue, military 
readiness, and personnel training. Two basic categoriss of 


requirements originate fron these responsibilities. The 
first type will be termed SHIP-SPECIFIC requirements. These 


assignments arise from requisites physically associated with 
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a particular ship and her crew. There is a need to maintain 
her equipment and machinery, and train her crew to function 
as a unit so that the ship can perform general operational 
missions. For example, Refresher training (Reftra), Navy 
ASW exercises (Navex), and Maintenance periods (MAINT) are 
SHIP-SPECIFIC requirements. Each ship must be scheduled to 
individually satisfy these needs. The second type of 
requirements, termed GENERAL, cover the direct operational 
areas of responsibility. Successions of ships are used to 
collectively satisfy these requirements. Alaska fishery 
patrols (Alpat), Academy Cadet training, and oceanographic 
data ccllection (Ocean) are exanples or GENERAL 


requirements. 


The geographical area where the mission is perforned is 
also an important consideration. Travel time to and from the 
mission area 15 necessarily present in the one-day 
resolution of the published schedule. (For instance, travel 
time from Hawaii to Alaska is 7 days.) The time spent by a 
ship away from its homeport is affected by this travel time. 
Also, in order to minimize travel time to the mission area, 
different missions are performed by the Hawaii-based ships 
than by the continental-based (CONUS) ships. 


E GUIDELINES 


The GENERAL and SHIP-SPECIFIC requirements are specified 
in three ways by the scheduling guidelines. The first is 
the setting of the number of ships that are to be on-scene 
in a given area at a given time (e.g. During May, two ships 
Should be near the Aleutian Islands). This type of 


Specification is called a PATROL requirement. The second is 


the setting of a repetitive frequency to be met (e.g. once 
per quarter). The fulfillment of this type of requirement 
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usually takes a standard amount of time and is to Бе 
accomplished without interruption- The third is the 
specification of a total amount of time to be spent in a 
given period. This time quota will usually be divided into 
several segments (e.g. The maintenance requirement of 13 
weeks per year may be divided into segments or 4, 4, and 5 


weeks each.) 


Three principal  morale-related guidelines affect the 
amount of time each ship can be scheduled to be away from 
its  homeport. First, a maximum limit is set for the total 
"Away Homeport" (AHP) time per year for each ship. Second, 
it is desirable to balance AHP time between ships since it 
Mnc cions as a pseudo-measure of each ship's share of the 
workload. Third, a maximua limit is set for the duration of 
any single cruise (tne time from departure from honeport to 


Boturn). 


Other guidelines concern the iongest and most difficult 
GENERAL requirement - Alaska Fishery Patrol  (Alpat), A 
maximum limit on the total yearly Alaska Patrol time рег 
Ship is set. If a ship is assigned back-to-back Alpat 
missions, the minimum iintervenung in-homeport pericd is 
eight weeks. Before any Alpat, a minimum four week inport 
period is desired for adequate preparation. It is also 
desired to alternate snips that are assigned Alpat missions 


in the rough weather months of winter. 
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III. CHARACTERIZATION OF PROBLE! 


ШО CHARACTERIZATION 


Examination of the published schedules for the period 
January 1975 to December 1976 snows the following basic 
elements: 6 ships, 2 year time span, 6 areas of 
responsibility generating requirements, and 129 sevarate 


missions used to satisfy the requirements. 


The frequency of mission-to-mission transitions made by 


the ships is shown in Figure 1. 


ee ee ip ee сс me END un ee A A E O MM 


975 to December 1976 


=e Sk ee cp -> = a ee E ъз шт = Se енә 


СОЕ БОВЕ ЕВЕ лое 
1 


шароле Alpat Oceani Reftra Naver Maint 


From: 

Aport 3 12 7 5 D 
Alpat 11 4 1 
Ocean + 14 

Reftra 7 1 
Navex 6 1 1 
Maint 15 2 2 y 1 

Figure 1 
FOL instance, ships nade the transition from 


Oceanographic data collection (Ocean) to Alaska Patrols 
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(Alpat) fourteen times in the two years. 


The schema used by Conway, Maxwell, and Miller [3] to 


describe scheduling problens is a useful classification and 


helps 


to succinctly define which factors are known and 


unknown by the scheduler: 


1. 


or 


Requirement Arrival Process - At the start of the 
Scheduling process, the total requirements for each 
area of responsibility are specified within a given 
expected range. The number of separate missions that 
will be necessarv to fulfill each requirement is not 
Specified DUE OE to the Scheduling process. 
Requirements are not simultaneously available but 
become available at individual times according to 


frequencies and timings given in the guidelines. 


Resources Available - The number and individual 


— m coms co Sc 


capabilities of all ships are known by the scheduler. 


BNcWwEPattern - As shown іп figure 1, the mission 
transitions that occur have definite patterns. The 
Sequencing of missions is critical, mostly because of 
tne limitation on cruise length. Also, the PATROL 
requirement introduces precedence and linking 
relationships between the individual nissions used to 
fulfill the requirement: the end of a patrol for one 
are sieuid implicitly specify the start of a patrol 


Er another ship. 


Measure 


O «р 


previously, there is no single MOE but rather a complex 


(us Effectiveness (HOE) - As explained 


combination of requirements and goals. 


This problem also possesses constraints beyond the scope 


those given by Conway, Maxwell, and Miller. These 


constraints are caused by the physical limitations of tne 
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Ships and the needs of their crews. The total AHP time 
scheduled for each ship should be equally balanced with 
respect to the other ships. The length of single cruises 
should be less than the maximum limit or the endurance of 
the ship or crew. For each ship, the total Alaska Patrol 
time per year should not exceed the specified limit. These 
three constraints, referred to as the "limiting" constraints 
ШІ Тигеһег discussions, are not typical production or job 


shop constraints. 


To summarize the Coast Guard problem, the scheduler has 
definite resource information and definite knowledge cr the 
type of transitions that can occur. He knows the total 
mission area requirements and must determine the number and 
durations of the separate missions that collectively satisfy 


these requirements. The requirements are not simultansously 


available at the beginning of the period but arrive 
Бааша са11у. The "limiting" constraints impose additional 


restictions on the sequence and durations of separate 


Missions. 


This proplem is very large. The final schedule has 
one-day resolution for about 5,000 ship-days and about 6 
Dission-types. Any economically feasible analytic model 
nust aggregate and divide the scnedule period into weeks or 


Even months and then schedule in terms of these units. 


DAS ELECTION OF APPROACH 


An Integer Linear Programming approach has been 
considered. This approach to scheduling problems (as 
demonstrated by Prabhaker [10]) 1s not feasible because of 
the large scale of this Coast Guard problem. À schedule for 


6 ships, 6 mission-types and 2 year time span with 1 week 
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Fesoluticn would yield about 3,600 0/1 variables. This size 
is too large for commercially available computer codes. A 
special code for this particular application would have to 
be develcped. (There is no apriori guarantee that even the 
best possible approach would solve tne problem with a 


reasonable amount of computer resources.) 


Simulation has been discarded as a solution method for 
Several reasons. Long development time is required to 
implement the proper sequence of random starting solutions 
and priority or heuristic-guided search metnods. (Panwalkar 
and Iskander [9] list over 100 static and dynamic scheduling 
ENedispatching rules.) This approach is typically expensive 
to run. Difficulties exist in interpreting the resuits: data 
sensitivity and experimental variation are difficult to 
differentiate. Effective external controls for guidance or 
model coersion are usually lacking. There are also serious 
reservations about the potential Flex iD day of the 
resultant computer program to meet future problen changes 


(in the guidelines, for example). 


А foundation of optimization theory and techniques is 
desired. The general approach initially selected to address 
this Coast Guard scheduling problem is the hybrid Quadratic 
Assignment/Linear Programming model of Geoffrion and Graves 
[4]. The detailed description of this model is presented in 


the next Chapter but a brief summary is appropriate here. 


A  Quadratic Assignment model determines a candidate 
Sequence or separate missions for each ship. A Linear 
Prcgramming model then determines the duration of each 
mission. Then, two heuristic searches vary the candidate 
sequence and new durations are obtained in an attempt to 


improve tae schedule. 
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This model has been selected for the reasons listed 


below: 


1. 


The flexibility in the problem statement accomodates 


nicely the indefinite information on the number and 


Sduration Of individual missions. 


No simplifying assumptions or deletion of guideline 
restrictions appear to be necessary, allowing an 


uncompromising view of the entire scheduling problem. 


The Quadratic Assignment model naturally nandles the 
sequencing restrictions (which normally cause great 
difficulty and computational complexity for analytic 


models). 


The continuous Linear Programming model narurally 
hongles the "irmiting" constraints of this problem and 
allows the lengths of the missions to Бе optimally 


adjusted to pest meet the guidelines. 


The Linear Programming model ( uSing flexible penalty 
functions) is a natural way to incorporate the 
guidelines and  morale-related goals into a linear 


Segective runction. 


The model handles dynamic, non-simultaneous arrival of 


requirements. 


The model pcssesses special structure amenabl=2 to 
implementation techniques appropriate for the large 


scale of this problen. 


This hybrid model has been successfully implemented Гог 


es industrial preduetion application of smaller scals. 
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OVERALL PHILOSOPHY 


The model proposed by Geoffrion and Graves makes several 
significant contributions to solving an interesting a class 
of scheduling problems: the separation ОТ the 
coabinatorially dum сЕ Sequencing problem fron the 
continuous ailccation and timing problem; a flexible problem 
statement reflecting true managerial discretion; and a 
method to express this class of scheduling problems in the 
База, finite formulation required for the Quadratic 


Assignment problem. 


The class of scheduling problems addressed by this 
lud method is stated in production terminology as 


fcllows: 


There are several "similar" continuous process 
facilities (lines) operating in parallel. Each is 
able to manufacture some subset of products with 
пажа production rates and costs. Significant 
costs are incurred by each changeover of a line 
мәс рссацсага One product to another. Production 
orders are received dynamically and are nor 
sinultaneously ауа lable for scheduling., Each 
ОЕ ШЕ Ка срес су с tOtcimanoOunt of production 
to occur between an early start time and a late 
finish time. An order can be split among lines ог 


produced non-contiguously on the same line. The 
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desired schedule has the minimum total production 


and transition costs over a specified horizon. 


The solution approach followed by their mcdel allows the 
broadening of this problem statement to allow lower and 
upper limits on the order demands rather than an exact 
amount. ICO улак оп Of the start and finish times of 
orders, and early or late completion by individual lines is 


allowed but penalty costs are incurred. 


In the hybrid model, a Quadratic Assignment (QA) problem 
is formulated and solved using the exact demands, the 
Specified horizon, and the early start and late finish 
memes, Violations of the horizon or start/finish times incur 
penalty costs. The QA solution is not an optimal, minimal 
East solution. The problem is solved by heuristic methods 
Mat obtain very good, low cost final results that are 
"locally optimal"! (i.e. certain easy changes made to the 
solution cause increased cost). Рог the remainder о? this 
thesis, "minimal" implies this locally optimal condition. 
The resultant assignment (lines to orders over the fixed 
time span) is then used only for the sequence of crders 
assigned to each Production line; Each continuous 
Production run on a line, called a “campaign,” 91s noted. 


The quantity and timing information is discarded. 


Given the product sequence rrom the QA solution, a 
Linear Program (LP) is generated and solved to find the 
duration of each campaign on each line, using the bcunded 
demands, the desired early start and late finish times, and 
a flexible schedule horizon. The LP minimizes total 
production and penalty costs (incurred by violations of the 
desired order times and/or horizon). The sequence costs are 
determined by the QA sequence and are not changeable by the 
EP. Next, two local searches are employed to examine the 


sequence of campaigns and locate potentially better product 
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Sequences. For each favorable candidate sequence, another 
ИР solution is obtained. The final solution is that 
combination of sequence and product durations with the 


minimum total production, sequence and penalty costs. 


Be QUADRATIC ASSIGNMENT FORMULATION 


Although many variants of the Quadratic Assignment 
problem have been presented in the literature, this 
discussion covers the models influencing the Coast Guard ОА 
model. Discussed are the original Koopmans-Beckmann model 
[7], the generalization addressed by Graves-Whinston [6], 
the Specialization Of Graves-Whinston used ру 
Geoffrion-Graves in their scheduling model, and (in the next 
Chapter) tne large scale QA model arising in the Coast Guard 
prcblem. The reader will note that the specializations made 
in the Coast Guard model are appropriate to this particular 
proplen. The techniques to handle the large scale can be 
applied to generalizations of the Coast Guara model with 


some additional effort. 


The original statement of a Quadratic Assignment 


(QA) problem was made py Koopmans and Beckmann for the 


assignment сі n "indivisible" manufacturing plants ton 
fixed geographical locations 50 that inter-plant 
transportation costs of опе commodity are minimized. As 


Stated by Lawler ( 8], let 


€ - transportation cost per unit from 


оса оп | го Location k; 
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Ч. = quantity shipped from 
a 
E piantao plane p (interaction quantity); 
mimo assigned tO location j, 
È 
x 0 otherwise. 
The object is to ninimize the interplant 


meansportation costs, 


n n n 
E x x : 1) 
> > Mus | ў кі mr ' 
р= 1 i=] j 


Eurject to 


n 
X... 

і = 1 

n 
х = 1 EOD i 
1) 
= 
Es, the optimal pairs of (plant i, location j) are 


desired. 


Graves and 4hinston have generalized the 
Koopmans-Beckmann model by the addition of a fixed cost tern 
Bene objective function. The assignment of plant 1 to 
ШОО Ооп 1 сап incur Known, fixed costs (for example, 
purchase of land, installation of highways and other 
services). The problem is to minimize total transition plus 


E ced costs in determining the (plant, location) pairings. 





3. Geoffrion-Graves Formulation. 


The Geoffrion-Graves formulation is a specialization 
ЕЕ the Graves-Whinston QA problenm. Restrictions are 


МЕ ед on the "piant"-to-"plant" interactions that occur. 


In the G-G QA formulation, the scheduling horizon 
ООН) on each production line is divided into equal 
ndivisible time-units. Each time-unit is a "plant" in the 
Koopmans-Beckaann sense. The interaction between time-units 
(plants) affects only those time-units immediately 
preceeding and following them on the same production line. 


Mié” interaction О matrix (composed of all q ) has the 
ір 


structure Snown in Figure 2. 


Q-Matrix 
TIME-UNITS 


о 
о 
me Om e 


TIME-UNITS M5 


en айты ЖЕТ 
= 
о - 
о 


Й е 
> 


A en 
| 
TE 
| 
| 


Figure 2 


(Note: The Geoffrion-Graves model is generalized to handle 
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similar" BEoOdUctyon lines. ее са» product3on and 
transition costs between lines are not required. 
Proportionality between lines for these costs may exist. 
The +1 entries in the Q-matrix are replaced by the line's 


mmisiticn cost proportionality constant.) 


The quantity A of each production order is divided 
MO equal indivisible "product-units." Each product-unit is 


a Koopmans-Beckmann "location.". 


The finite guantization is determined by the choice 
of the basic time resolution S. Tne nuaber of time-units 


per line 1l is: 


3 
И 


Hu Ts (2) 


the schedule horizon. 


x, 
5 
(D 
ri 
(D 
EI 
li 


The number of product-units per production order k 


Шит а product p is: 


па в (кж су, (3) 
К К р 
where ze = tne order's specified demand 
R = the product's standard production rate. 
р 


Requiring that the total number of time-units equal 


the totai number of product-units yields: 
п = = m = > n. (4) 
I k 
| К 


The problem is infeasible if total product-units 
exceed total time-units. A "slack" product can be added to 


enforce (4). 
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The quadratic assignment problem can now be stated 
ШЕ) the set or pairings (time-unit i, product-unit j) such 
that total transition costs plus total fixed costs is 


minimized, 


n n n n n n 
а fC л хт а X, 
ip jx ij pk La 49] 
feel i = 1 EE ESSENT EI 


subject to (5) 
n 
x = 1 for j^ in 
ij 
1 = 1 
n 
х = 1 LOE тї = 1уп 
17 
хау 
where: 
С 2 transition cost between product-units j and к; 
jk 
q. = 1 if time-units i and p interact, 
1р 
= Q0 otherwise; 
d = fixed cost of assigning Product unit j to 
17 
time-unit i; 
ЖИЕ - Ut product-unit 7 Us assigned to time-unit i, 
1j 


0 otherwise. 


fbe Ci natrix, showm in Figure 3, Contains the 
EMSuuct-unit to  product-unit transition cost, derived from 
mee Standard product-to-product costs. Each product-unit is 
associated with an order which is for a particular product. 


(The standard transition cost between product p and product 
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p' is denoted by T 2) 


рр" 
спаз 
MISSION-UNITS 
n, | na 
| | | 
E | | 
n о ГТП | 
| | AP, AP, 
| BE ou 2 
| | | 
ДА | о | T | 
n | 
| | 
MISSION-UNITS ---4------------ 
| | | 
Т ЖТ о | 
БЕС Mors га ыға de 
| | | 
| T | | 
: | ! | 
| ¿ | | 
Figure 3 
Ec ОП пле D matriz өл Fixed costs. о Inis 


i) 
Esc contains costs (incurred when product-unit j is 


assigned to time-unit i) which account for the following 
efíects: 


IEEE The production costs; 


2. The transition cost, for each line, from the last 
product on the previous schedule to the first in this 


scheduling period; 


3. Prohibitive penalties for producing an order before the 
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Сл шас ыле or aftez che late finish time; 


ШИ Inrteasibr)irty of certain products which cannot be 
produced by certain lines (a high cost of infeasibility 


is used). 


(Note: There is an explanation of how this ОА 
formulaticn is interpreted for ships and missions of the 
Coast Guard problem in the next Chapter. Also, Appendix A 


contains an example and picture.) 


See GiNEAR PROGRAM FORMULATION IN G-G MODEL 


A change of variables occurs between the (time-unit, 
product-unit) pairs of the QA and the time durations cf the 
LP. Each continuous production run of a product 
(product-units for the same order) is noted and called a 
"campaign." The primary LP variables are the durations of 


each campaign. 


Without changing this sequence of campaigns on each 
line, the LP determines the duration (in continuous time) of 
each canpaign so that the minimal totai cost is obtained. 
EHE three cost components in the LP's objective function 


are: 
IE Production costs, 


IEEE enaltres for violating orders' early start or late 


finish time, and 


Bee Penalties for modifying the specified horizon, Н, оп а 


line. 


The total quantities for an order are constrained to fall 


between the soecified lower and upper limits. 
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№. IMPLEMENTATION. 


Geofírion and Graves have implemented their model to 
Schedule six chemical reactors for a one month planning 
horizon for Dart Industries, Inc. Various kinds of plastics 
are the products. The QA problem is solved by the Graves 
and Whinston method [ 6]. 


The Graves and Whinston method solves a generalization 
of the Koopmans-Becknann QA problem. То Simplify this 
discussicn (and without loss о? generality), the special 
Structure of the Q-matrix of the G-G model is assumed. The 
G-W method is an n-stage decision process where n is the 
total number of time-units for assignment. At the Rn stags, 
ШӘ! palrings of  (time-unlit +, roduct-unit j) have been 
made. In the rasic method, no backtracking occurs so these 

th 
X-1 assignments are permanent. Тһе К Stage examines all 
possible pairings of the n- (k-1) unassigned units. For each 
pairing, the expected final total cost, comprised of three 


components, is calculated. These three compcnents are: 
1. The costs incurred by the previous x-1 pairings, 


2. The immediate cost of assigning the Balz being 


th 
considered at this К stage, and 


3. The expected value of future costs that will be 
incurred by the remaining n-k unspecified assignments, 


given that the current pair is made permanent. 


The expected value of future costs is calculated 
explicitly py a formula rather than Бу enumeration. The 
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assumption is made that the assignment probability of each 
possible permutation of assignments is equal. The 
replacement of enumeration for unspecified assignment costs 
by the conditional expected value 1s a significant 
contribution cf the G-W method. 


The pair, designated by (1', j'), with the minimum 
expected total cost is selected by the assignment criteria. 
Mext, the fixed costs, d, , of those time-units that 

1) 
inmediately  preceed ог fellow i' are updated. For example, 
mete time-unit i preceed i'. Any product-unit j assigned to 
time-unit i will be followed by  product-unit j'. The 
EeEISIition cost of 7 to j' is now known and is thus a fixed 
cost. The fixed cost elements for  time-unit i are thus 
Wpdated by the addition of this known transition cost. A 


(ат update for the time-unit following i' is also done. 


th 
The methcd then proceeds to the k stage. 


Simce no conditional computations occur at each stage, 
the computational effort and memory requirements of this 
portion of the G-W method may be determined as a simple 


Bunctjion of n. 


А  lccal search (called Switch) is next performed cn the 
(time-unit, product-unit) pairings. All pairwise exchang2s 
of prodtect-units are cyclically tested "ior possible 
improvement in total cost. For each new pairing, the change 


Ene rixred costs and transition costs is deterministic. 


At the termination of this Switch algorithm, а ОА 
solution is obtained. Next, tne sequence of campaigns is 
extracted and made the incumbant sequence for the LP stage 
of the G-G model. A Linear Programming solution is obtained 


for this sequence. Then, two different local searches are 
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performed on the candidate sequence. The first search 
(called Slide) moves each campaign into all other possible 
positions. The second search (called Switch) makes a 
pairwise interchange of all campaigns. For each new, 
favorable sequence of campaigns selected by a screening 
criteria, an LP solution is obtained. If the total cost 
(sequence plus LP costs) is inproved, the candidate sequence 
replaces the incumbent and the two local searches are 
restarted. When both searches terminate, the incunbent 
sequence and durations constitute the "iocaliy" optimal, 


minal solution. 


The LP problems are solved by G-G with an LP subroutine 
described in Graves [5]. This algorithm includes advanced 
features for dealing with degeneracy, employs a complete 
mutual primal-dual simplex mechanisn, and uses a basis 
representation with a small, segregated inverse. The 
prcgram uses complete in-core storage of the coefficient 
matrix, and handles logical variables implicitly. Ranging of 


equations and bounded variables are not available. 
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Bee OVERALL CONSIDERATIONS 


Toe Coast Guard scheduling model is based on the 
Geoffrion and Graves approach. The Area Commander's 
Scheduling guidelines, expressing mission requirements and 
the morale-related goals and constraints, are converted into 
the objective functions and constraints or the Quadratic 
Assignment and Linear Programming models. The objective 
function expresses pseudo-costs that are incurred DY 
deviations from the reguirements and goals. The optimization 
models, by minimizing their objective functions, seek the 
closest compliance of the resulting schedule to the 
guidelines. The objective functions are summations, each 
term expressing with penalty costs the degree of compliance 


with one of the requirements or goals. 


In the Coast Guard problem, a ship corresponds to a G-G 
meme." The time available on a ship 15 divided into 
discrete time-units. Each mission requirement is a G-G 
"production order" and is similarily divided into 


nission-units. 


The guideline's requirements and goals have not been 
completely captured py the QA and LP models. Some reatures 
are included in the QA that are not in the LP and vice 
versa. The guidelines not modeled at all are left for human 


action upon the resultant schedule. 
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ЕРИ QUADRATIC ASSIGNHENT MODEL 


The modeling for the Quadratic Assignment algorithm is 
performed by setting the structures and components of tne 


three matrices: 


C - mission-unit to mission-unit transition 
costs, 
Q - interaction between time-units, and 


D - fixed costs. 


Each of these matrices are specially structured for this 
particular Coast Guard problem. The C-matrix will coerce the 
mission sequences to follow the actual patterns shown in 
Figure 1. The Q-natrix is the same as the G-G matrix. The 
D-matrix will express in fixed costs the desired 
start/finish times, the two types of mission requirements - 


PATROL and single continuous duration, and other itens. 


ab EE GS a eco E A oe 


The mission-unit to  mission-unit transition costs 
are directly determined from tne requirement-to-requirement 
transıticn costs. Frequently occurring mission transitions 
are strengly encouraged by low costs and undesired 
transitions are discouraged by high costs. (А transition 
between  mission-units of the same requirement incurs no 
cost.) Figure 1 in Chapter III shows the frequency of the 
various transitions in actual schedules. The mission Inport 
(ship in homeport on search and rescue standby) plays the 
dominant role. Most of the time, Inport is the mission 
occurring between all other missions. The exception, Ocean 
Шівсвапсагардіс data collection), із performed enroute or 


returning frcm Alaska Patrol. The number of occurrences of 


37 





each extended sequence is shown below. 


Inport/Ocean/Alpat 14 (approx) 
Alpat/Ocean/Inport 4 


AE ron ission A to Mission B" transition 
eoSt, a ПОЗЕ СН ЛЕУ arises in modeling these two 
alternatives. The goal of the model is to match the 
scheduler's logical ranking with the ranking determined by 
the pseudo-costs structure. The example below shows the 


problen. 


Eron-=io Wransıtjon’Dirfrieulty 


ча шек аа щш u сс a me En m m mn mn ne Damm 


Sample Transition Costs: 


Tor: 
Inport Alpat Ocean 
ЕШӘОЕТ 0 30 10 
From: Alpat 10 0 30 
Ocean 90 10 0 
Logical Pseudo-Cost Transition 
rank Rank Cost 
Sample Sequences с a ES 
Inport/Ocean/Alpat/Inport 1 1 30 
IUDOPt/Alrpat/Ocean/Inport 2 3 150 
iurort/Ocean/Iinport 3 2 100 


Figure 8 


The sequence Inport/Ocean/Inport would never be 
scheduled and thus is highly unfavorable. The pseudo-cost 
mechanism can not convey this as the difference in the 
rankings shows. The rrom A to B" Structure cannot handle 
"from à to B to C" sequences. A possible way to prevent the 


undesired and unrealistic Inport/Ocean/inport sequence fron 
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having an inappropriately low cost 15 to eliminate 
Ocean/Alpat or Alpat/Ocean. Based on the frequency Table, 


Ocean/Alpat is retained. 


The same difficulty occurs if Travel time is modeled 
as a mission to preceed or follow those missions requiring 
it. The irrational sequence Inport/Travel/Inport would nave 
a lower cost than a logically prefered sequence. (e.g. 
Inrort/Travel/Reftra/Travel/Inport). Thus, Travel time as a 
mission is not explicitly modeled in either the QA ог LP 
models. This exclusion is not considered serious because of 
the times (2 to 7 days) are constant and dependent upon the 
specific ship. Thus, allowance for tnis Travei time can be 


made in the duration of the Inport missions. 


The transitions that are assigned low (encouraging) 
Costs are listed. All other transitions are discouraged by 


high costs. 


Inport/All missions 
All missions/Inport 
Ocean/Alpat 


ces ann ARO A сощы: 


The Coast Guard Q-Matrix is the same as that of 5-6 
with the specialization that the transition cost 
Eroportionality constants equal опе for ali ships. This 
matrix serves to associate each time-unit with only the 
immediately adjacent time-units on the same ship (first and 
Hast time-units on a ship are special cases). This structure 
is the same as that of the Multiple Traveling Salesman 


prcblem in the open literature. 
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3. D - Matrix 


The D-Matrix contains the fixed costs of pairing 
time-unit i to mission-unit Jj. Five different fixed cost 
components are used to model five guidelines (The reader can 
refer to Appendix A where the initial D-Matrix for the 
Sample problem contains labeled examples of each of these 
components). The first three are the same as in the G-G 


model. 


The first component models the requirements that can 
be fulfilled by a specific ship. The geographical locations 
of homeport and mission area primarily determine possible 
asSignments for the GENERAL requirements. This component 
also models the SHIP-SPECIFIC requirements by allowing only 
the intended ship to be assigned its SHIP-SPECIFIC missions. 
If a ship cannot satisfy a requirement, the intersecticns of 
all associated time-units and mission-units are assigned 


very high costs. 


The second component provides the mission-to-mission 
transiticn ccst between the mission assigned to a ship 
immediately prior to the start of the current schedule 
period and the first mission it will be assigned in this 


period. 


The third component coerces the enforcement of early 
Start and late finish dates that are associated witn sach 
requirement. Each time-unit that follows the late finish 
time of a requirement is assessed a penalty cost for 
lateness. This ccst increases linearly as the lateness 
increases. Similarly, a penalty cost for preceeding the 
specified early start time is assessed. These assessments 


are nade EOE all mission-units associated with the 
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requirement. To model the predetermined nature of the first 
three to six months of the schedule, and the freedom of the 
last six months, the rate of linear increase of the penalty 
can be contsobled over the time horizon for each 
requirement. A requirement in the beginning of the schedule 
can be assigned a higher penalty rate than a requirement at 
the end. Thus, a deviation of the same amount can be 
penalized differently, depending upon the proximity to the 
beginning of the schedule and the flexibility allowed for 


the requirement. 


Ше толЕсп  соцропепе of the D=Matrix does not 
correspond to one in the G-G model. Some requirements 
Stipulate that the desired duration of the reguirement 
cannot be split and must be completed at one time. This 


Е т= usually a SHIP-SPECIFIC requirement. As the QA 


nodel allows splitting of requirements into multiple 
missions, this fourth component attempts to prevent 
splitting by imposing a Structure on the pertinent 


(time-unit, mission-unit) pairs. Figure 5 shows with and 
without States for a small example. Ап >additicnal 
interpretaticn is placed on the mission-units; they are to 
be assigned ina left to right order (the left-most first, 
and the right-most last). As only one nission-unit can be 
assigned to each time-unit, the first time-unit after the 
early start time of the requirement is encouraged to have 
the leftmost (first) mission-unit assigned to it. The 
time-unit preceeding tne late finish time is also encouraged 
to have the rightmost (last) mission-unit assigned to it. 
Similar encouragements are used as appropriate for the 
remaining mission-units. The encouragement is done by 
appropriate displacement of the time penalties of the third 


component. 


The effect of this imposed structure is to decrease 


the number of pairings that incur no time penalty cost. The 
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(time-unit, sission-unit) pairings are thus encouraged to 
follow one of the "paths" shown by the arrows. This 
structure does not guarantee that the requirement will not 


be split, but does penalize a split. 


Single Mission Structure 
Without With 
Mission-Units Mission-Units 
1 2 3 1 2 3 
1 1 1 1 1 2 
2 Bar lyiStart 
Time- 3 
Units 4 1 
5 Late Finish 2 1 
6 1 1 1 3 2 1 


(Note: The numbers are units of Time Penalty.) 


Figure 5 


ИСТ ГЛ component Of the D-Matrix explicitly 
models the PATROL requirements (given number of ships on 
scene at same time). А sSinilar left to right time 
orientation is placed on the mission-units. By the nature 
of this type of requirement, a strong relationship is 
thereby created between particular time-units. That is, the 
first mission-unit should only be assigned to the one 
time-unit of each ship that corresponds to the same actual 
calendar time. Figure 6 depicts the following situation: 
for a schedule horizon of three weeks (with one week 
resolution), two ships must satisfy a PATROL requirement of 
1 ship always on scene. Thus, the second mission-unit of 
the requirement is strongly associated with only the one 
time-unit of each ship that corresponds to the second week 
of the period. All other time-units are discouraged from 


possible assignment by high costs (denoted by M). 
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4. Guidelines Not Incorporated in QA 


Those parts of the guidelines specifying the 
mOrale-related and "limiting" constraints are not modeled in 


the QA model. No mechanism has been found to condition tae 
e zum en 
assignnent of the k stage on, for example, the amount of 


AHE time previously assigned each ship. These constraints 


are incorporated in the subseguent Linear Programming model. 


C. LINEAR PROGRAMMING MODEL 


The Linear Programming nodel for the Coast Guard 
: У \ z ^ч "73 
problem contains the same oasic structure as tne G-G model, 


Additional penalties are incorporated in the objective 
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monet lon ЕС modei the "limiting" and morale-related 
guidelines. The problem statement of the QA model is 
broadened. The single target duration of a requirement is 
extended to have lower/upper limits and the time horizon can 


vary on each ship. 


The primary variables of the LP represent the 
mission durations. (Note: For the remainder of this thesis, 
a mission refers to a G-G campaign.) All other variables 
are non-negative penalty variables that measure the 
deviations from desired conditions. The six types of penalty 


variables in the objective function are listed below. 


1. For each ship, the difference between the sum of all 
Missions on the ship and the time horizon minus the 


starting time of the ship. 


2. For each ship, the difference between the sum of all 
Away Homeport missions and the specified total desired 


amount for that ship. 


3. For each ship and each mission-type with a maximum 
limit, the positive difference between the sum cf all 
missions on the ship for this mission-type and the 
Specified limit. 

4. For each missicn, the negative difference between the 
starting time of the mission and the associated 
requirement's early start time, and the positive 


difference between the actual start time and the late 


start time specified. 


5. For each mission, the positive difference between the 
start time plus the duration of the mission and the 


associated requirement's late finish tine. 


6. For each sequence of Away Homeport missions, the 
positive differnce between the sum of durations in the 


cruise and the limit on the duration of a single 
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роте. 


The constraints of the model, with no provision for 
violation at any penalty cost, pertain to the upperylower 
limits on the total requirements and the durations of each 


mission. 


1. For each requirement, the sum of durations of the 
missions associated with the requirement is within the 


lower and upper ranges. 


2. For each mission, its duration is between specified 


lower and upper bounds. 


The PATROL requirements are not handled explicitly. 
This type of requirement is handled implicitly by the 
BE LOLIty assigned to the  start/finish times of the 
requirement. The priority weighting of the linear penalty 
rate is the same as that explained in the  Quadratic 
Assignment model. Guidelines such as the minimum of eight 
weeks inport between successive Alpat missions and the 
minimum cf four weeks prior to an Alpat for preparation are 
explicitly handled Ьу manipuiation of the mission duration 


bounds. 


DESEGUIDELINES NOT HODELED 


Some facets of the Coast Guard problem have not been 
analytically modeled. Travel time is not explicitly handled 
because of the sequencing difficulty already discussed 
previously in this chapter. The target total Away Homeport 
time of each ship could be adjusted to partially account Гог 
the effect of travel time. The single cruise limit could 


also be adjusted. 
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Multi-ship conflicts are not modeled. Overlar of 
first/last weeks. or Reftra, overloading of districts with 
ships simultaneously an maintenance, staggering the 
Hawaii-based ships, and alternating winter Alvats are 
exanples. These must be checked and possibly corrected by 


the scheduler. 


The desire to have at least one ship on Inport (search 
and rescue standby) at all times has not been modeled. 
Presently, the ship workload 15 such that this 15 
practically guaranteed. The effects of the new 200 mile 
ИШПИ ПО lirit may prevent this Inport coverage from 
occurring naturally. Inport can then be handled as a PATROL 
reguirement missicn. District training and non-scheduled 
operations are not modeled since they are not manually 


scheduled at present. 


Single day timing of homeport, return, departure, and 
on-scene dates must be left to the manual scheduler. 
Overloading of the fuel pier in Kodiak, Alaska, and Sunday 
arrivals at Reftra are examples of considerations that the 
schedulers handle when producing the single day resolution 
of the published schedule. 


The items not analytically modeled can be included in an 
edit program that would examine the resultant schedule for 
these relationships and inform the schedulers of hand 
erections to be nade. Such a utility program could also 
prove very useful for computing "pseudo-costs" of alternate 


schedules manually produced. 
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vi. IMPLEMENTATION 


AA A SÍ we TS A AS A 07 ыш, | шр 


Pee OVERALL CCNSIDERATIONS 


The G-G model and computer program described in [4] 
provide an outstanding example of successful implementation 
pa c complex optimization model for production use. The 
model concisely captures the pertinent features of the 
prcduction scheduling problem in an optimization context. 
The computer system efficiently provides good solutions to 
the model for the intended client. Unfortunately, the sheer 
size of the Coast Guard problem makes the G-G system 
hopelessly expensive to use. Accordingly, a new system has 
been designed for the large scale Coast Guard problem which 
preserves all the exciting model features of G-G, while 
incorporating additional model features and yielding vastly 


inproved computational efficiency. 


The primary considerations for the implementation of the 
analytical model are: to overcome the formidable computer 
memory and execution time difficulties arising from the 
large scale of the problem, and to test and evaluate the 
proposed analytic model. The Geoffrion and Graves code 
(FORTRAN IV) is used as a starting point. It is designed 
for commercial use for smaller scale industrial production 


scheduling problems as discussed in Chapter IV. 
An ISM 360/67 at the W. R. Church Computer Center of the 


rs SP Naval Postgraduate School has been used for 


development work in this study. 
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The overall structure of the implementation of the Coast 


Guard analytical model is depicted in Figure 7. 


B. QUADRATIC ASSIGNMENT IMPLEMENTATION 


The Quadratic Assignment code used for an initial start 
is an igplementation of the algorithm developed by Graves 
and Wninston. The code uses explicit  in-core storage for 
cach of the three matrices: р, Es ando. Thus, data 
storage for a 1,000 unit problem would be approximately 3 
million wcrds which is clearly an uncomfortable size. Thus, 
the first obstacle to overcome to enable the solving of a 


1,000 unit quadratic assignaent problem is core size. 


The G-W algorithm is implemented for the most general 
case of the QA problem. With the general Q-matrix, many 
conplex relationships between "plants" can be specified (the 
preceeding/fcllowing relationship of time-units is a 
specialization). The work factor (approximation for numbers 
of computer operations or expense, as a function of probien 
size) for the general algorithm is estimated to be on the 
order of 8. By specialization of the algorithm to the 
defined structures or the C and Q matrices of the Coast 
Guard anad tocco model, reduction of core and improved 


performance will result. 


The (multiple traveling salesman) specialization used 
for the О matrix is to collapse the n*n matrix into two 
vectors cf size n. One vector acts as an index set to give 
for each row the column number of the single non-zero 
element (or Zero if none is present). The second vector 


performs the inverse mapping from each column to the гом 
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with the non-zero element, if any. Figure 8 depicts these 
Ewcc vectors for the Q-matrix illustrated in Figure 2 
(Chapter IV). 


Row i: 1 2 3 ... nh -2 m -1 р 
Cclumn # cf Non-Zero: 2 3 4 To n -1 fa 0 
1 1 
Necrosgtor Columns ot Q-Matrix LforgEirst Ship 
Column i: 1 2 3 sss тел! п 
1 1 1 
Row # of Non-Zero: 0 1 2 "TE Dres n di 


Figure 8 


For instance, for row 3, the non-zero element appears in 
column 4; for column 3, the non-zero element appears in row 


DE ROW ШЕ has no non-zero element as it is the last 
time-unit on the ship (no time-unit follows it). Column 1 


has no non-zero element (the first time-unit on tne ship has 


no predecessor). 


2 
This change reduces storege requirements from n to 2n, 
3 


ü 
EUM he work factor from n ton. 


The special block structure of the C matrix enables the 


elimination of this matrix resulting in a storage reduction 
2 . ‘ . 
of n. By using several levels of indexing from the 


nission-tc-mission transition costs, the appropriate 


mission-unit to mission-unit transition cost is obtainable 
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from data arrays already used elsewhere in the model. 


The final storage reduction is obtained from the 
observation that at most two updates are made to each 
element of the D matrix. After a (time-unit,mission-unit) 
assignment is made, the fixed cost D matrix is updated. The 


mieurted transition cost into or out of the assigned unit is 


added to each d, іп the rows preceeding and following the 
1) 


row just assigned. Thus, the maximum value a d, element 


Can attain is the maximum initial fixed cost plus twice the 
maximum transition cost. Ву scaling the initial fixed costs 


aad һе transition costs such that the largest possible d... 
1) 


15 
is less than 2 (32768 is absolute magnitude).  Integer*2 


(16 bit, nalf-word) storage is used rather than integer*4 


Dt, rfull-word). Integer arithmetic operations are 
about 10% slower than floating point arithmetic on the 
machine used in this study, but the data storage reduction 


has to be made. 


By these changes, the total data storage requirement for 


2 
the three matrices is reduced from ЗЕ words to 
| 2 у 
ЕО: васе1у п /2 words. Thus, a problem with 1,000 time 


units would require 500K words rather than 3,000K words for 


data storage. 


From trial runs, three additional algorithnic 


specializations have been noticed for tne Coast Guard 
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Beoplen. Mission-units for Inport comprise опе fourth to 
one third of the total number of  missicn-units. All 
mission-units of Inport have identical initial fixed costs. 
All dynamic transition cost updates maintain this identity. 
Thus, during the expected value calculation matching each 
unassigned tine-unit with all unassigned nission-units, only 
one of the unassigned Inport mission-units need be tried as 


a candidate rather than all n/3. 


The second specialization arises from the observation 
that alacst all of the Inport mission-units are always the 
last mission-units assigned to time-units and that the 
expected value of the future cost of thesé new assignments 
is zero. Once the point is reached where only  Inport 
nission-units remain unassigned, the identical nature of all 
of these units means that the expected future cost of their 
assignment becomes permutation independent. These units can 


be assigned to any unassigned time-units and in any order. 


The third specialization depends on the likelihood that 
infeasible (strongly undesired or impossible) assignments 
occur. Ыс to i calcula ing the expected value of 


assignment of (time-unit i, mission-unit j), the fixed cost 
for this assignment, d, , can be tested. If this fixed cost 
i 


has exceeded a predetermined infeasibility value, the trial 
assignment can be preemptively skipped. A difficulty arises 
if an infeasible assignment should really take place via the 
lowest Dean criterion. With this efficiency-motivated 
modification, abnormal termination  witn an incomplete 
assignment map would then occur since the algorithm would 


bypass tnis infeasible assignment. 
These three specializations have not been made for 


reasons apparent later in this Chapter (a new algorithm for 


the QA model is developed). It is speculated that their 
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affect will reduce the work factor of the specialized G-wW 


3 3 
method frem n ton /4 if Inport missions comprise one third 
of the schedule. 


Following the completion of a (time-unit,mission-unit) 
map by minimum expected values, an additional algorithn 
(called "Switch") is tried to obtain quickly and easily 
further improvement in the objective function. All pairwise 
interchanges of the mission units are cyclically tested by 
exhaustive enumeration. For this Coast Guard problem, 
striking improvements are achieved. This effect and other 
basic perforrance data of the G-W method is listed in 
Figure 9. (The 5 day time resolution has not been run 


because the execution time would be excessive.) 
C А NEW QUADRATIC ASSIGNMENT METHOD 


The Graves-Whinston method, with the modifications and 
Specializaticns described above, is estimated to use four to 
eight hcurs of computer time (IBM 360/67) to solve a 1,000 
unit Quadratic Assignment problem. Since this is excessive 
(especially considering that the purpose of the overall 
analytical mcdel is ultimately to aid manual schedules in an 
interactive fashion), a second method of implementation has 
been explored. Other reasons for introducing a second 


method are: 


IENEThe ")imiting" and morale-related constraints are not 


modeled in the QA; 


2. The QA sequences wiil be subject to increasing 
modification as the LP model incorporates these new 


types of constraints; 


3. Significant cost reductions are made by the Switch 
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transformation; and 


4. The fixed cost matrix for the Coast Guard problem is 


heghly structured. 


For these reasons, the new method is designed to quickly 
obtain a feasible assignment of (time-unit 1, mission-unit 
j) pairs using only the fixed cost D-matrix. А Linear 
Assignment problem is solved to obtain a set of  pairings 
With the minimal total fixed cost. Then, the fixed cost 
matrix is upcated with transition costs in the same manner 
as the G-W method. The same Switch algorithm is then used 
Bosminimize total fixed and transition costs. The final 
solution, obtained at the end of this algorithm, is a 


Poemueron tO the original Quadratic Assignment problem. 


The Linear Assignment problem is solved by a special 
subroutine called CGNET, which is an adoption for tne Coast 
Guard problem of the well-known primal simplex network 
package SNET [1]. The package is modified to use the 
completely dense assignment (fixed) cost matrix and to 
exploit where possible the special structure of these 
extremely large problems. Several external parametric 
controls permit tuning of the package for efficient 


performance. 


The network problem for the 5 day resolution has 860 
nodes and 739,600 arcs. CGNET constructs an optimal Linear 
Assignment solution in 11.3 minutes (ІВМ 350/67). A 
literature Search indicates that this is the largest 
assignment problem for which an optimal solution has been 
eonstructed. The previous record was a 450,000 arc model 
done in a feasibility study for the Navy Computer Assisted 


Distribution and Assignment (CADA) model. 
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The memory requirement of а. words for data storage is 
the same as the G-W method's. Figure 9 shows the 
comparison between the G-W and CGNET methods for computation 
time, effectiveness and computational cost of the Switch 
algorithm, and total pseudo-cost of the final (time-unit, 
MS on- unit) pairing. Note that CGNET produces excellent 


quality sclutions with significantly less computational 


expense. 

Performance of Quadratic Assignment Methods 
Resolution 20 DAY SIDA 10 DAY 5 DAT 
Size (n) 206 288 420 860 
Method е! CGNEI G- W CONST G-W CGNET G-W ЕСИЕТ 
Assignment 
Time (min) 85225 02380020278 (jo) 221510, 243.7 == TE: 
Switch 
Time (min) 0227 0.42 0752 0.92 18955 1.91 -- 7.66 
Total est 
Tine (min) 8.50 0.80 21.30 1 5298 70555 4.28 480. 18.98 
Assignment 
Cost 81110 148730 96120 172740 110060 182070 -- 168230 
Final Cost 37470 52570 38590 36480 45960 40940 -- 11100 
€ Switch 
Improved 53.8 64.7 59.9 78.9 58.2 19,5 -- 93.4 

Figure 9 
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ШИ ЕОЗТ-РВЕОСЕ55ІЗС6 OF THE QA SOLUTION 


The Graves-Whinston and CGNET methods for the QA problem 
do not guarantee that the optimal solution will be found. 
In addition, the solution ma y contain infeasible 
assignments, multiple campaigns when only one is desired, 
and illegal transitions. To determine if these undesired 
events have cccurred, quick and effective edit checks are 


conducted. 


in the hybrid approach of the analytic model , the 
Quadratic Assignment solution is examined and used £or only 
the pure sequence cf missions conducted by each ship. The 
consecutive time-units of a ship are checked for each 
sequence of misSion-units associated with the same 
requirement. Checks are made ror assignment of missions to 
ships for requirements that they cannot fulfill. Any such 
mission that is found is relocated to be the first mission 
on the first ship able to conduct tne mission. The 
structure placed in the D-matrix to encourage a single 
Mission for a requirement when this is desired does not 
guarantee this result. Thus, an edit check that retains the 
first missicn encountered and discards the rest is 
performed. Another quality not guaranteed by the QA 
solution is the absence of all infeasible or undesired 
transitions, such as Alpat/Reftra. An optional edit check 
can eliminate this type of transition by the insertion of an 
OT t mission wherever needed. However, as later 
procedures may cause the reintroduction of such transitions, 
this check is best delayed until after the final LP solution 


r obtained. 
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E. IMPLEMENTATION OF LINEAR PROGRAMS 


In this hybrid model, once the mission sequence is 
obtained, a linear program is solved to deternine the 
mission durations that minimize the objective function. 
Then, neighboring mission sequences are generated and tested 
to see if improvement occurs. A neighboring sequence is 
defined as any Sequence that is generated by applying one of 


the following operations to the present sequence: 
1. Move one mission to another position (called Slide). 


2. Interchange the position of two missions (called 


Switch) . 


This trial of neighboring sequences is absolutely critical 
for this Ccast Guard problem because of the differences in 
the content of the QA and LP models for this problem. The 
QA soluticn's sequence will be far from optimal because the 
"Limiting" ccnstraints are not considered. Many neighboring 
sequences will be favorable candidates to reduce total cost 
and for each such neighboring sequence, an LP solution is 
necessary. Thus, the computational reguirements, core and 


time, must be minimized for each LP call. 


The initial LP code (used in the G-G implementation) 
Solved a problem cf 330 rows and 330 variables in the range 
Of 17 to 21 seconās on the IBM 360/67. Analysis of the 
constraint equations shows that only 1% of the total 
coefficient matrix is non-zero and the only non-zero value 
Ше ғ1. Thus, an LP code using a data structure for the 
Sparce constraint coefficients is dzsired. Also, а code 
allowing a reduction in the number of equations by ranging 


and bounding methods reduces the dimensionality of the 
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problem and thus reduces computational tine. 


Since the QA and LP portions of the hybrid model are 
distinct, the program and data areas of each model can be 
overlayed in main memory, resulting in significant savings 
of core. With a projected LP problem size of 300 rows and 
500 variables, the original LP code would have storage 
requirements larger than either specialized QA code for time 


resoluticns cf 20 days or less. 


A new LP system called XS has been installed for use 
with the Ccast Guard problem. XS [2] is a prototype 
optimization system designed to serve as a research testbed 
for evaluaticn of advanced design features for large scale 
linear, nonlinear and integer problems. Тһе system is used 
as a subroutine, with the LP problem set up provided by the 
calling discipline. Salient features of the version of XS 
employed include implicit ranging (upper and lower values 
for each ccnstraint), logical bounding (upper and lower 
values for each variable), and effective external controls 
for tuning the system to perform well оп the class of 


problem at hand. 


The LP cbjective function, described in Chapter V, is 
composed cf  non-negative penalty variables whose values 
measure the deviation of the schedule from the desired 
goals. The deviation is obtained by expressing the goal as a 
constraint, with the insertion of the appropriate slack or 


surplus penalty variable in the equation or inequality. 


28 





For exanple, the deviation of a ship from the desired 
total AHP time is determined by adding the following 


eonastraint tc the LP: 


Objective Function: (cost * slack) + (cost * surplus) 
Constraint: 

Total All 

AHP = AHP t Slack -~ Surplus 

Е Missions 


Thus, the izplementation of the LP model has tne two 
types of non-violable constraints and six penalty-related 
constraints (that determine the value of the penalty 
variables). Several cf these types of constraints аге 
independent of the overall number and structure of the 
sequence cf missions. These are the total horizon time, AHP 
time and Alpat time on each ship, and the early/late start 
Ete finish tine penalty constraints. The number of 
time penalty equations is fixed by the number of missions 
for requirements possessing these penalties. Since no 
missions are deleted or added by the LP, no change in the 


number of missions occurs. 


The numrer of equations for the cruise duration penalty 
and the quantity bounds on the requirements can уе reduced 
for sach individual LP based on its particular mission 
sequence. The number of Away Homeport cruises is dependent 
on the individual sequence. For each sequence of two or 
nore AHP missions, an equation is introduced. Rather than 
introduce an equation for a solitary AHP mission, the 
nission duration's upper bound is compared to the AHP cruise 
Et апа lowered to this value if necessary. This method 
does not allow a solitary mission to exceed the cruise 


limit, even for a penalty assessment. (This is a change to 
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the strict definition of the LP model given in Chapter V.) 


The last constraint type is the quantity restriction on 
each requirement. These are hard constraints with no 
provision fcr penalties for non-compliance. An equation is 
introduced only if a requirement is split into two or more 
missions. If only one mission produces the requirement, all 
previous bounds on the mission duration are supercedsd by 
the minimum and maximum quantities of the requirement. Note 
that if the single mission forms a solitary AHP cruise and 
the maximum quantity exceeds the AHP sequence limit, no 
penalty will be charged if the cruise limit is exceeded as 


no equation for cruise penalty is generated. 


The typical LP fora 6 ship, 2 year schedule has 171 
rows and 435 variables. The number of rows for 2ach type of 
constraint is: 6 each for total horizon, total AHP time and 
moral ALPAT time; 86 for early/late start times; 51 for late 
finish times; 12 for cruises of 2 or more missions; and 4 
requirement duration totals. There are 145 mission 


variables and 290 penalty variables. 


The nature and balance of the objective function is 
Critical to obtaining the desired type of schedule. The 
balance between the six types of penalties will determine 
the relative importance and preferences of the tradeoffs 


that will be made during the minimization process. 
The cost structure for each type of penalty is also an 


important factor. Figure 10 depicts the present 


structures. 
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Е. IMPLEMENTATION ALTERNATIVE POR LP 


Examination of the solutions to all test LP's shows that 
mission durations were aiways exact integers. This leads to 
the investigation of a Network formulation оЁ the LP. A 
complete transformation to a network has not been found for 
the CG model because of a few complicating constraint types. 
The horizon limit, time penalty and requirement quantity 
constraint types (which collectively comprise the entire 
original G-G model and about 85% of the equations of this 
Coast Guard model) can be rerormulated as a pure network 
nodel. This reformulation has not been implemented, but a 
significant reduction in time for each LP call, especially 


for the G-G model, would occur. 


The variables (Arcs) for the network formulation are the 
length and the starting time for each mission. Figure 11 
pictorally shows the fundamental structure of the network. 
There is cne equation 9 (node) гог each mission of a 
requirement that has any of these qualities: a late 
completicn penalty, a following mission with an early start 
or finish date. Each requirement with multiple missions 
will have an equation. Single mission requirements are 
handled ry appropriate bounds on the mission's length. The 
total time limit equation for each ship is added between the 


last missicn node for each ship and the root node. 


The tEalancing of AHP time is handled in the netwcrk by 
its complement, the time spent on Inport missions. Three 
types of missions emerge from the root node. The first type 
1s all AHP time requirements. The second type is all 
SHIP-SPECIFIC homeport requirements. The third type is all 


GENERAL ship in-homeport requirements, usually only tne 
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Inport (SAR standby) requirement that is the system's slack. 
To obtain the total non-AHP time on a ship, the type two and 
three missions on a ship need to be added in some manner. 
If the type two missions do not have a fixed total, this can 
not Бе done ky the network. If they do have a fixed limit, 
appropriate limits and penalties can be placed on the sun of 
type three nissions going to each ship. Actually, 
calculations to check for this situation шау not be 
worthwhiie in actual implementation. Replacement by one 
non-network equation per ship is simpler. Non-network 
equations are also necessary for each AHP cruise sequence of 
2 cr more missions and one per ship for each mission-type 


with maximum limits (i.e. Alpat). 


An important property of the network is тһе resulting 
pictorial display of the model. Based on this picture, a 
method to perform dynamic changes to the formulation may be 
developed tc eliminate the complete problem generation for 
each new LP mission sequence. Alternatively, a mission that 
16. а candidate for relocation may be placed in several 
places at once. Sensitivity analysis and reoptimizaticn can 
determine the one best placement faster than a succession of 
complete LP calls. This and an ability to reoptimize based 
on the current solution may drastically reduce the time and 


increase the information gained from гасп LP call. 
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б. CALLING STRATEGY FOR THE LINEAR PROGRAMS 


The first sequence of missions obtained from the 
Quadratic Assignment solution andthe post-QA/pre-LP edit 
processing will be far from optimal because the "limiting" 
constraints are first introduced in the LP model. Thus the 
ability to improve the solution by evaluating many 
"neighboring" mission sequences is vital to the successful 


sclution of the scheduling problem. 


The first algorithm (called Slide) performed on the 
mission sequence removes eacn mission in turn and inserts 1+ 
into all other possible locations. The number of dirferent 

2 
mission sequences to be evaluated in some manner is n. Ап 
actual call for a full LP solution, that is, complete 
explicit enumeration, is clearly inpossible. д сана 
strategy implicitly evaluates each possible sequence and an 
LP is solved only for the most promising candidate 
sequences. Additionally, after a neighboring sequence is 


found with an improved total sequence and LP penalty costs, 
* > * * 2 
it becomes the incumbent solution and another n cases ale 
evaluated until no further improvements can be found with 
this transform. 

The second algorithm (called Switch) evaluates pairwise 
interchanges cf missions. The number of possible switches 

. . 2 

for each incumbent sequence is (n -n)/2. Again, after an 


inproving switch is found, the transform is restarted. 


For each type of transform, a calling strategy to decide 
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ШШШ ТО invest in an LP solution is critical. Those 
transformed sequences that introduce assignments on ships 
incapable of fulfilling the requirement can easily be 
detected. But for feasible assignments, the expected change 
in each of the six types of cost penalties must be 
calculated or estimated. The change in the sequence 
transition cost is deterministic and easily calculated. The 
changes in the penalties, however, have an unpredictable 
nature because they are dependent on the mission lengths 
which are set by the LP to obtain the minimal objective 


Function. 


To solve this problem for Slide, a detailed calling 
strategy is implemented for five of the penalty types. (The 
horizon penalty is not included.) The transition cost is 
deterministically found and estimates are made fcr the 
change in start and finish time, total AHP, total Alpat, and 
АНЕ cruise penalties. Amazing accuracy has been achieved. 
In the early stages of the algorithm, a high percentage of 
the LP calls tnat are made lead to improvements extremely 


close to the prediction. 


Other methods to maximize the contribution of each LP 
Call are tteoretically similar to pivot DIEI mg and 
Selection rules. For Slide, a mission is moved to each 
feasible new position and the expected cost change is 
calculated. The proposed new position with the most 
beneficial expected change is tne only one considered 
further. The predicted change is then compared to a decision 
rule cutoff value. Based on this test, an LP call is either 
made or the Slide algorithm discards that mission as a 


candidate and proceeds for the next mission. 
In another strategy, the decision cutoff value is 


changed as the algorithm progresses. Presently, the sliding 


decision rule has four values, each higher in value (-8000, 
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-4000, -1000 and -600). Those sequence changes that cause 
very large and beneficial cost improvements are actively 
sought in the opening stages. After one decision rule value 
leads to a complete circuit of Slide with no improvements 
being found, the rule assumes the next higher value and the 
algorithm is restarted. The computational time spent doing 
these calculations and tests for a complete circuit of the 
algorithm is less than half the time of a single LP call. 
Figure 12 shows the dramatic improvements made in the 
Objective function, the percentage of LP calls being made 
that lead to true improvements, and the overall success of. 


these calling strategies for Slide. 


For the second algorithm, Switch, a prediction of the 
cost change is made for only the start and finish time 
penalties and the total AHP time penalties. The transition 
cost change is deterministic and, thus, is calculated. Only 
one value of the decision rule cutoff value is used. LP 
solutions are obtained for all favorable candidate sequences 
rather than for only the most favorable. The performance of 


B Salgorithm is also indicated in Figure 12. 
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Performance 
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ОЕ OF THE MODEL 


Before discussion cf the qualitative properties of the 
schedules produced py the implementation of the Coast Guard 
analytical scheduling model, the general types of input data 
and user external controls must be understood. There are 


seven types cf input data used by the model: 


1. Number of ships, number of mission types, the time span 
of the schedule, and the time resolution of the ОА 


model. 


2. Description of the mission types: away/in nomeport, 


single (no split of requirement), and/or Patrol. 


3. Description of the ships: capability for each mission 
type, last mission on previous schedule, initiai fixed 


portion of present schedule. | 


4. Description of requirements: mission type, target 
demand for QA formulation and lower/upper limits for LP 
formulation, early and late start times and late finish 
ше priority of tines (fOr control of tine penalty 


rates). 


5. Pseudo-Costs: mission-to-mission transition costs and 
cost rates for the penalty variables (Horizon, 
start/finish times, total AHP time, total Alpat tinme 


апа cruise duration). 
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6. Desired goals for total AHP time, total Alpat time, and 


TIMISec duraticrn. 


7. Implementation controls: LP calling strategy decision 


EHt-Off values, level of output, length of maximum run. 


By the specification of the input data, the scheduler 
has effective control of the resulting schedule. The cost 
rates for the penalty variables determine the balances and 
tradeoffs that will be made by the optimization process to 
arrive at the final proposed schedule. By adjusting a rats, 
the influence and importance of that guideline's goal is 
changed with respect to all others. An example of cost 
balancing between transition costs and the cruise  penal-y 


rate is given in Figure 13. 


Another important control is provided by the cost 
penalty structures shown in Figure 10 (Chapter VI). The 
influence this structure has on the tradeofís and 
preferences within each guideline goal is illustrated in 
Figure 14. A structure with larger penalties for large 
deviations forces the AHP balance closer to the desired 


goal. 


A powerful property of the hybrid QA/LP approach is the 
separation of the combinatorial sequencing of. missions for 
each ship and the determination of the duration of each 
mission. This approach and the implementation techniques 
provide the capability of bypassing the QA, using a 
previously generated sequence, setting the penalty rates and 
obtaining a potential scnedule. (Note: Since the QA model 
does not include the un ng and morale-related 
constraints, changes іп the penalty rates or desired goals 
for these constraints affect only the LP portion of the 
nodel.) Dynamic interactive use of the model is an inherent 


property cf the modeling approach. 
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Two Alternatives for Distributing Total AHP Time 


Sr: Ship 1: four units under desired balancing amount 
ӘНЕС 2: опе over 

Б: Ship 1: two under 
Shop 2: three over 


6 4PENALTY 
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deslred. total amount. 


Figure 13 
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Situation: 


Alternatives: 


А 
В 


ca GNNED-WERe cup Er mm 


Between Two Components 


Schedule 10 weeks Alpat 

weeks Ocean | 
Ignore all penalties except Cruise 
limit and Transition costs. 
Cruise iimit is 10 weeks. 


Alpat and Ocean are AHP missions. 


Inport/Ocean/Alpat/Inport 
Inpcrt/Ocean/Inport/Aipat/Inport 


Cost Structure I 
Transition Costs 
inport Alpat Ocean 
Inport 0 10 10 Cruise penalty 
Alpat 10 9 1000 2000/unit over 10 
Ocean 1000 10 0 
Ccsts: Transition + Cruise = Total 
A: 30 + 4000 = 4030 
B: 10 30 + 0 = 1030 
B is preferred but as discussed in a previous Chapter, 
Inport/Ocean/Inport is an infeasible sequence. The penalty 
Cost of Ocean/Inport should be more than that necessary to 


link Ocean/Alpat even i£ cruise limit is violated. 


cost Structure II 
Transition Costs 
inpcrt Alpat Ocean 
Inport 0 10 10 Cruise penalty 
Alpat 10 0 1000 2000/unit over 10 
Ocean 7000 20 0 
Costs: Transition + Cruise = Total 
: 30 + 4000 = 4030 
В: 7030 Б 0 = 7030 
А is preferred, as desired by the guidelines. The cost 
of Ocean/Alpat was raised above the increase їп cruise 


рше у caused by adjoining 2 


pat. 


weeks of Ocean and 10 weeks of 


Figure 14 
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The model can be used for both strategic апа tactical 
decisions. Strategic "what 1?" questions can be examined 
for their effects on scheduling (for example, changes in the 
number of ships available, the cruise limit for a new class 
of ship, or new areas of responsibility for the Coast 
Guard). Dynaric schedule changes, due to requirement changes 
or unanticipated ship non-availability, can be evaluated and 


the least "upsetting" solution found. 


PROUAL ITY OF SCHEDULES PRODUCED 


Because cf the importance of managerial judgement and 
personal preferences, no single MOE is available to 
quantitatively measure the quality of the schedules prcduced 
by the analytic model. 


The Schedules produced conpletely satisfy the 
restricticns governing desired and undesired mission 
transitions. No transitions with high, infeasibl2 costs are 
present. Almost all requirements are scheduled to meet 
their specified start and finish times. Those that are not 
Alpat missicns are usually one or two weeks displaced from 
the desired times. The Alpat missions, having the Patrol 
requirenent, are exceptions. As explained, the Patrol 
requirement is implicitly modeled by setting high priority 
(causing a higher penalty rate) on fulfilling the specified 
start/finish times. About 70% of the Alpat missions are 
Scheduled tc satisfy the times, but the remaining Alpat 
missions vary from being displaced 1 week to a maximum of 8 


weeks. 


This result for a Patrol requirement mission-type is not 


realistic. To oktain even this result, it is necessary to 
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express the cverall requirement as many detailed, individual 
requirements with exact durations (equal iower/upper bounds) 
and specific start/finish times. This level of input detail 
is an undesirable feature. This result and the necessary 
detail of data are implementation issues, and not model 
shcrtcomings. model. An enhancement discussed in Chapter 


VIII will correct this poor performance. 


With the exception of the Patrol type of mission, the 
present model produces good quality solutions. The 
guidelines that are explicitly modeled and those that can be 
implicitly modeled by the scheduler through the input data 
capture all of the major considerations of the overall 
guidelines. This model also meets the objective of quickly 
generating several potential scnedules. Minor variations in 
some control parameters, particularly the minimum and 
Maximum durations for Inport missions, can cause alteration 
of the schedules. For the Coast Guard problem, this ability 
to easily generate alternatives is a highly beneficial 


property. 


«Тһе model, in summary, provides good quality solutions, 
generates alternative schedules easily, and is highly 
responsive tc managerial controls. This model can be a 
powerful managerial tool for scheduling the Hign Endurance 


Cutters. 


The use of this model for Medium Endurance Cutter 
scheduling (not the direct focus of this study) appears more 
in doubt: principally because shorter ship endurances cause 
shorter patrol lengths; more multi-ship relationships exist 
for Search and Rescue coverage; Reserve training cruises are 
highly complicated; and larger relative lengths of travel 


time to patrcl time exist. 


74 





ПОО / LP RESULTS 


Ише Size and computational difficulty of the Quadratic 
Assignment problem is determined by the time resolution. 
The time resolution (S) is the unit that determines the 
number of finite time-units to be scheduled for each ship, 
and the number of mission-units for each requirement. The 
results displayed in Figure 9 (Chapter VI) show the 
computaticnal performances of the Graves-Whinston and CGNET 


methods rcr time resolutions of 20, 15, 10 and 5 days. 


Compariscn of these methods and/or resolutions cannot be 
based completely on the final ОА solution cost. Between 
resoluticns, the number of  mission-unit to mission-unit 
transitions is different, causing differences in transition 
cost. Between methods, the final cost is an indication of 
the degree of compliance of the schedule to the desired 
mission sequences and start/finish times. TIMES HOGAR 
indication of the quality of the overall solution that will 
result at the end of the LP since the "limiting" constraints 
are not modeled in the QA. Also, the contribution of the QA 
1s only the mission sequence extracted Irom its (time-unit, 


wsio unit) pairs, not the additional timing information. 


The CGNET method usually produces a larger nunker or 
Separate missions than the G-W method. Most ог the 
additional missions are Inport which is beneficial later in 
the "neighboring" sequence searches. Having more Іпрог%5 


causes larger transition cost reductions to be possible. 
Qualitative evaluation of tne model's final schedule 


shcws that the method of solution of the QA has no impact on 


the final quality. Each method produces a superior final 
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БЕШЕСІ Ле fOr about half of all cases.  CGNET, with its good 


computaticnal performance, is the recommended method. 


Contrary to expectations, the 20 day resolution prcvides 
"higher" quality overall schedules than the smaller (5 day, 
10 day) resolutions; the 5 day resoluticn produces the worst 
Schedules. Тһе schedules from the finer resolutions have 
greater difficulty in satisfying the transition restrictions 
and reducing the violations of the desired start/finish 


times. 


The computation time for each LP solution is primarily 
determined Бу the number of missions. Aith about 145 
missions, an LP requires, on the average, 0.5 seconds for 
problem generation, 4.5 seconds for solution, and 100 
Pivots. Тһе LP calling strategy implemented is of vital 
importance - the most significant improvements are located 
and made first. In this way, the total number of individual 
LP solutions 1S minimized. In all experiments, the 
schedules with the lowest pseudo-cost require fewer total 
calls to the LP, and have the highest hit ratios. This 
result means that after the first 20 to 30 significant 
improvements, extended computer time and more LP calls are 


not particularly effective. 


To summarize the computational performance, the best 
schedules are obtained by the CGNET method for tne QA with 
20 day time resolution. Computation time for this method is 
about 52 seconds. With various bounds on Inport missions, 
the LP portion of the model requires an average of 7 
Minutes, 36 seconds. Thus, good potential schedules are 
obtained by this analytic model in an average 8.5 minutes 


using 70 K words of memory. 
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A. APPROACH 


With tne overall success of the hybrid approach, major 
changes in the modeling approach are not  necessarjy. With 
the exception of the PATROL requirement, the items discussed 
in this Chapter are ideas for possible improvement of 
computaticnal performance and additional managerial ccntrols 


that would be useful in a production level implementation. 


Improvement (fewer total LP calls) would result if the 
"limiting" and morale-related constraints can be considered 
in the QA. The inclusion of these constraints requires a 
nechanism to condition each assignuent (on the total number 
of AHP  time-units previously assigned each ship, for 
example). The Graves-Whinston method, which uses 
probability in determining the next assignment, can 
potentially te extended. So far, the necessary theoretical 


developments in this area have not been made. 


The Linear Prograuming aodei needs explicit modeling of 
PATROL requirement missions. А печ type of penalty or 
non-violabie constraint linking the end of one ship's patrol 
and the beginning of the next ship's should be added. (The 
next section gives further discussion of this point.) Also, 
additional equations for AHP balancing for each year as well 
as the whcle schedule horizon may be needed. Avxdairftaculc 
decision is where to split the mission sequence into tnose 


AHP missicns of the first year and those of the second year 
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prior to the mission durations being known. 


The intrcducticn of human interaction at several points 
nay alsc be worthwhile. Visual inspection of the sequence 
from the QA solution, as well as periodic inspection of the 
pregress of the LP calling strategy and sequence searches 
шау improve computation times and solution quality. The 
БЕСИ оп Of a post processor to inform the scheduler of the 
status of those guidelines not analytically modeled would 


also be appropriate for a production level progran. 


B. IMPLEMENTATION 


Several implementation modifications will increase 
Managerial ccntrol and further enhance computation speed. A 
more complicated selection scheme for deciding which one of 
Multiple missions to retain for requirements that cannot be 
Split would remove at least 10 to 15 LP calis. Also, 
relocaticn of  infeasibly assigned missions to the best 
rather than the first location would be an inprovenent. А5 
Shown Бу Figure 14 (Chapter VII), a change to a more convex 
penalty structure for total AHP time would lead то tetter 
workload balancing between the ships. Travel tiae can be 
explicitly modeled in the LP model. The travel times 
between areas for each mission transition are known. The 
lower/upper bounds on mission durations and the desired 
total time for each ship can be dynamically calculated 
during the generation of each LP at no increase in LP 
Solution time. All of the above listed ideas are easy 


modifications to the present implementation. 
A more difficult implementation modification is the 


explicit modeling of missions with the PATROL requirement in 


the LP model. The alternate network formulation discussed 
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in Chapter VI most easily accomodates this change. Іп this 
formulation, each mission's starting time and duration are 
available as primary variables. The linking of one 
mission's finish to another mission's start is easily and 
directly acccmplished. This enhancement will remove the 
only difficulty with the quality of the final prcposed 
schedules, and also remove the undesired detail presently 
required in the input specification of PATROL requirement 


nissions. 
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APPENDIX А 


SAMELE PROBLEM AND RESULTANT SCHEDULE 


Problem: Two ships are available for a 7 week schedule 


horizon tc fulfill these requirements: 
1. Seven weeks of Alpat with PATROL requirement; 
2. Two weeks Maint for Ship Two between weeks 5 and 6; 
3. Two weeks of Ocean between weeks 2 and 5; and 


4. Three weeks of Inport. 


Alpat and Ocean are Away Homeport missions. Only Ship 
One can fulfill the Ocean requirement in this period. The 
last mission on the previous schedule for Ship One is 
Inport; Ship Two, Ocean. The time resolution for the ОА 
model is one week. The following list contains the cost 


penalties and desired goals: 


1. Infeasible assignment cost 700. 
2. Time penalty cost rate 50% 
3. AHP penalty ccst rate 600. 
4. Cruise penalty cost rate 2000; 
5. AHP desired goal рег ship 4 weeks 
6. Cruise limit 5 weeks 


Transitticn Costs: 


INPORT ALPAT OCEAN MAINT 
IBNPORT 0 10 10 10 
ALP AT 10 0 700 700 
OCEAN 700 10 0 700 
MAINT 10 700 700 0 
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mni tial D-Matrix of Fixed Costs 
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